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Based on the work of Chen, Lii and Pope, we derive expressions for the _D > 6 dimensional 
metric for Kerr-(A)dS black holes with two independent rotation parameters and all others set 
equal to zero: ai 7^ 0, 02 7^ 0, as = 04 = • ■ • = 0. The Klein-Gordon equation is then explicitly 
separated on this background. For _D > 6 this separation results in a radial equation coupled to two 
^ generalized spheroidal angular equations. We then develop a full numerical approach that utilizes 

I the Asymptotic Iteration Method (AIM) to find radial Quasi-Normal Modes (QNMs) of doubly 

rotating flat Myers-Perry black holes for slow rotations. We also develop perturbative expansions 
for the angular quantum numbers in powers of the rotation parameters up to second order. 
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I 1 I. INTRODUCTION 

I After the advent of the brane world scenario [1] and the AdS/CFT correspondence [2], there has been a growing 

interest in the study of higher-dimensional black holes. Perhaps the strongest driver behind this interest is that the 
pH Large Hadron Collider (LHC) [3] may produce black holes in the near future and such black holes, if produced, would 
I— I have large angular momentum. Therefore higher dimensional generalisations of the Kerr solution are the natural 
setting for these studies. 

Rotating black holes in higher dimensions were first discussed in the seminal paper by Myers and Perry [4] . One of 
^ I the unexpected results to come from this work was that some families of solutions were shown to have event horizons 
QQ for arbitrarily large values of their rotation parameters. The stability of such black holes is certainly in question [5, 6], 
p\J but no direct proof of instability has been provided. Another new feature of the Myers-Perry (MP) solutions was that 
T-H they in general have [ "^~"'" J spin parameters, making them somewhat more complex than the four dimensional Kerr 
solution. The first asymptotically non-flat five-dimensional MP metric was given in [7]. Subsequent generalizations 
to arbitrary dimensions was done in [8], and finally the most general Kerr-(A)dS-NUT metric was found by Chen, Lii 
I and Pope [9]. 

In the literature the focus has been largely directed toward solutions with only one rotation parameter, the so called 
^ simply-rotating case. The reason for this is that in the brane world, collider produced black holes would initially only 
• ^ have one dominant angular momentum direction. This is due to the particles that produce the black-hole being 
/\ confined to our 3-brane and therefore the rotation would be largest in the plane of collision on the brane. However, 
^ there is reason to believe this picture may be too naive. In any realistic situation the brane would be expected to 
have a thickness of the inverse Plank scale. At impact, the colliding particles would in general be offset in these thick 
directions and therefore further non-zero angular momenta would be present in other rotation planes of the black 
hole. Even though these angular momenta would be small compared with the rotation on the brane there is strong 
evidence [10] to suggest that such black holes would evolve into a final state in which all the angular momentum 
parameters were of the same order. There are also compelling theoretical reasons why one would want to go beyond 
the simply-rotating case. In particular, the Quasi-Normal Modes (QNMs) of these solutions may have applications in 
the AdS/CFT correspondence. 
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The study of the wave equations in higher dimensional rotating black hole spacetimes was initiated in [11] by 
analyzing the Klein-Gordon equation in five dimensions. The analysis relied crucially on the method of the separation 
of variables. The problem of separability of these wave equations in higher dimensions is a difficult one, even for the 
Klein-Gordon case; early attempts were only aimed at special cases [12]. Finally, using the Chcn-Lu-Pope metric, 
Frolov, Krtous and Kubiznak [13] were able to separate the geodesic equation and the Klein-Gordon equation in the 
most general setting. This was then realized to be due to the presence of hidden symmetries, in the form of Killing 
tensors [14-16]. A whole tower of Killing tensors and symmetry operators [17] can be constructed with the help of 
the corresponding Killing- Yano and conformal Killing- Yano tensors. They guarantee the separability of the geodesic 
equation, the Klein-Gordon equation and also the Dirac equation [18]. Unfortunately, the higher spin wave equations, 
especially the gravitational perturbation equation, have not yet been subjected to such an analysis. 

Even for the scalar wave equation, efforts so far have been focused mostly on the simply-rotating case. Notably, in 
[19-21], the stability of the scalar perturbation in six and higher dimensions was considered in the ultra-rotating cases, 
where no instability was found. Due to the interest in AdS spacetimes, the investigation was extended to Kerr-AdS 
black holes [22-26]. Here the expected superradiant instability did indeed show up. In addition. Hawking radiation 
in these spacetimes (Kerr-dS) [27, 28] were calculated with possible applications to the production and decay of LHC 
black holes [29]. 

In this article we show how both the D > 6 two-rotation metric and the separation of the scalar wave-equation on 
this metric can be achieved quite economically by working with the general Kerr-NUT-AdS metric described in [9]. 
That is, we begin with the full set of ["^^^J rotation parameters, Oq, and [ "^~^ J coordinate variables, yi, and then 
take the limit that all but two of the rotations goes to zero. This then reduces the full metric down to that with 
only two non-zero rotation parameters and allows us to present this metric explicitly in these coordinates, with an 
expression valid for any dimension D >6 ^. 

In this form the separation of the Klein-Gordon equation proceeds analogously to the case with all rotations present. 
For simplicity we will set the NUT charges Lq = although in general this is not an obstacle to separation. 

The structure of the paper is as follows: In the next section (Section II), we present the general metric of Kerr-(A)dS 
black holes with two rotations. The corresponding Klein-Gordon equation is separated into one radial equation and 
two angular ones. In section HI we develop a full numerical method using the Asymptotic Iteration Method (AIM) to 
solve for the angular eigenvalues and the radial quasinormal modes (QNMs) for two rotation parameters. Conclusions 
and discussions are then given in Section IV. In Appendix A the small 01,02 expansion for the angular quantum 
numbers for two angular equations analytically are given up to second order; we call this "double perturbation 
theory" . This allows us to check the AIM against the perturbative method. 

II. THE METRIC AND THE SEPARATED EQUATIONS OF THE KLEIN-GORDON EQUATION 

WITH TWO ROTATIONS 

In this section, we shall first derive the Kerr-(A)dS metric with two rotations from the general metric obtained 
in [9]. To specialize to the case with only two rotations, we take the limit 03, 04, • • • — > 0, while, without loss of 
generality, keep oi > 02 > 03 > 04 > • • • , where the a^'s are the rotation parameters. Then we explicitly separate the 
Klein-Gordon equation on this metric. The separation in the general Kerr-(A)dS metric was first performed in [13]. 
We follow their procedure for the case with only two rotations, where we find one resulting radial equation and two 
angular ones. 



We start looking at the metric with D = 2n, that is, for even dimensions. However, the result we obtain at the 
end is also valid for odd dimensions. For D = 2n, there are at most n — 1 rotation directions so we have a^, with 
i = 1, 2, . . . , n — 1. This metric, satisfying R^^ = —Zg^g^^, can be expressed as follows [9]. 



The case with D = 5 is exceptional because even though there are two rotation parameters there is only one Jacobi coordinate variable 
j/i . As we shall see, for D > 6, in the limit of keeping only two rotation parameters a total of two Jacobi coordinates survive. The 
D = 5 case will be considered in a separate work [30]. 



A. 



General metric with two rotations 
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where 

n— 1 n— 1 

U=Y[(^' + yl), U^ = -{r'+yl) n iVp-yl)^ (2.2) 

n—1 n— 1 

^ - (l + 5V)n('^' + «fe)-2Mr, ^c< = -(l-5'ya)n(«fe-2/')' (2.3) 

n—1 n—1 

w ^ n (1 - ff'y' ) ' = n («' - ) ' (2-4) 

a—l a—1 

71—1 n—1 

t = iHil-g^a^) , ^, = a,{l-g^a|)4>^ H {a^ - 4) , (2.5) 

i—1 k—l,k^i 

with 1 < a, z < 71 — 1. is the azimuthal angle for each a^^. With the direction cosines /x^, z = 1, . . . , ri, the metric 
for a unit _D — 2 sphere is just 

n n— 1 

df^^^^d^i^^^A^^d^^^ (2.6) 
subject to the constraint J27=i ~ ^- This constraint can be solved in terms of the unconstrained latitude variables 



UlJiial-yl) ^2 ^ UlJi yl 

a? Illlik^^{al - ai) ' " nr=i' 



Mi — 9T-rn-l / 9 9\ ' — T-rn-1 9 ■ 



a. 



Expressed in terms of the unconstrained coordinates the metric for the unit sphere becomes diagonal: 

n—1 n—1 



a=l 

with 



This choice then allows for a more symmetric form of the general Kerr-(A)dS metrics [9] as we shall see below for the 
two rotations case. 

To obtain a general metric with two rotations we take the limit 03, 04, a„_i — !■ while assuming that 
ai > 02 > as > • ■ ■ > a-n-i- From the definition in Eq. (2.7) we see that yi is of the same order of magnitude as a.;. 
Bearing this in mind, in the above limit, we can consider the different terms in equation (2.1): 

"Idr- = ill±MM±yll,r^ , (2.10) 

X Ar 

where 

A. = (1 + g'r'){r' + al){r^ + a^) - ^ . (2.11) 

Note the solutions of — lead to the black hole and cosmological horizons: r_|_ and r^^ respectively, for the Kerr-dS 
case (g < 0). 



^ Note that the metric for the odd case is slightly different, see [9]. Chen et al. also define an extra parameter an = in the even case in 
order to make some parts of the treatment between even and odd cases homogeneous. However, for better clarity we have elected not 
to do this here, i.e., we assume there are only n—1 parameters, ai, in the even case. 
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As another example, under this hmiting procedure, the quantities 
both seem to vanish in the hmiting process, however, the ratio 



C^3 , 2 
-^3 



\a\al) 



1 



yl 



dy. 



( yjyl 



gady. 



(2.12) 



(2.13) 



is actuahy finite. Here in defining we have taken into account only the angular variables associated with 03, 04, 
. . . , ttn-i- In the same way, part of the sum in the second term of the metric in Eq. (2.1) then constitutes 



n-l 

H^dyl 

»- n, 



a=3 



Similarly, part of the other sum in Eq. (2.1) gives 



E 

a=3 



(l + ffV)Ty ^- 
1 - q'^vJ, 



E 



f ylvl 

\alal 



af-yl 



E 9°'dyl 



a=3 



H 2 



Combining these two we obtain the metric for a _D — 6 sphere, 



/ 2 2\ /'n—l n—1 
^ 1 ^ / \a=3 i=3 



(yM] 



dQ 



D-6 > 



as indicated in Eq. (2.8). 

Finally, for dimensions D > 6, the metric with two rotations can be given by 



{r^ + yt)ir^ + yi) 



(1 - g'yf)il - g^y|) ^ - ytm ' yl) #i 

(l-g2a2)(l-g2„2) (l-g2„2)(„2_„2) 

{al-y\){al-yl) dc/i 



A 



yi 



ir^ + ytm-yi) 



(l-.g^ai)(ai-ai) a2 

(1 + ,9^2) (1 - g2y2) ^ (^2 + ^2)(„2 _ ^2) 



(l-.92a2)(i_32„2) (i_^g2„2)(„2_^2) „^ 



(1 - .9^ai)(ai - a?) ^2 



+ 



(r2 + y|)(y2_y2) 



dt 



{l-g^al){l-g^aj) {1 - g^aDiaj ^ 4) ai 



(r2+a2)(,2„^2) 
2„2^ 



(l-3^ai)(ai-a?) 02 



(r2 + y2)(r2^y2)^^^ , ('^' + y?)(?yi-2/?) 



2, ,2 ^ 



(yM) 



A,, 



dyj 



(r' + yDiyf-yl) 



A 



dy2 



V2 



and 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



Ar = [1 + g^r'^){r'^ + al){r'^ +al)-2Mr 
Ay, = {l-g'yl)ial~yl){al~yf) , 



7-D 



A,. 



{i~gV2){<-yi){4-yi) 



(2.18) 
(2.19) 
(2.20) 



Here we have kept the variables yi and t/2 instead of writing them in terms of angular variables. This is because the 
relationship, as shown in Eq. (2.7), is rather complicated to write out explicitly. If we solve t/i and 2/2 in terms of /zi 
and /Z2, we have 



2/1,2 = 



1 



2(1 - /i?) + ^2(1 - ^,1) ± J4a2a2(;.2 + ^2 _ 1) + („2(i _ ^2) + ^2(1 _ ^2))2 



(2.21) 
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It then follows that yi and 1/2 must be constrained by 

aa < 2/1 < ai ; < 2/2 < 02 (2.22) 

in order for Eq. (2.7) to be well-defined. 

We used a similar procedure as above to show that the metric in the odd dimensional D = 2n + 1 case reduces to 
the same form as the one shown in equation (2.17). 



B. Separated equations for the Klein-Gordon equation 

The separation of the Klein-Gordon equation in the general Kerr-(A)dS metric has been achieved in [13]. Here we 
shall show explicitly how the separation goes for the case with two rotations. To begin with it is convenient to rewrite 
the metric in Eq. (2.17) as 



-Qi 



Af^d^o + A\''d4'i + A\'>di^2 



,(1) 



(2), 



+Q3 A^°'^d^|Jo + A\^'>d^i + A]^>d^: 



(1). 



if' 



+ Q2 
2 



A^2^d^po + A'i^dijji + yl^^'dV'2 



Qi 



-dr^ 



Q: 



-dyl 



-dy. 



dnl_^ , (2.23) 



where 



Qi 



and V'fc are related io (pi and (j)2 by 

i>2 



Q2 



or conversely, 



a? 01 



Q3 



{r^ + yl){yl-yl) 



al(j>2 



g't 



{l~g^al){al-al) {I - g^al){al - al) 

ai4>l , 0202 



{l-g^aj)il-g^al) {1 - g^afKaj - aj) {1 - g^a^,){al ~ aj) 
gH (f>i 02 



{I - g^aj){l - g'^al) ai{l - g^al){al - a^) a2{l ~ g^aDiaj - af) 



ri 

02 

a2 



tpa + {of + aljijji + a^a20'2 
g'^tpo + (1 + g'^aDi-'i + a2ip2 



g20o + (l + .g al)V^i +alV': 



(2.24) 

(2.25) 
(2.26) 
(2.27) 

(2.28) 
(2.29) 

(2.30) 



The matrix A^'^'' is given by 



1 y1 + y2 ylyl 



1 

1 -r^ 



+ yl -r'^yl 



yl -r^yf 



(2.31) 



with the inverse B^^^y 



-Vl 



-2/2 



(fc) 



{r-^+vl){r^+vT) {r^+yi){yl-yi) {r^+yi){yi~yi) 
922 

Vl Va 

(r'^+yDir-'+yi) {r''+yi){yi-yi) {r'^+yi){yj-yi) 



V (r^+yt){r^+yi) {r^+yt){yi-yt) {r^+yi){'y{-yi) j 
In this notation the inverse metric components are 



(2.32) 
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9 



yiyi 



9 



V2 y2 



Q3 



9 



-^4)^0) 



^ r2 r)2 



Q: 



'^{i)^(j) ' (2.33) 
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plus the angular part related to the metric gab, and the inverse g , for a unit {D — 6)-dimensional sphere S . The 
determinant of the metric is then given by 



detg^^ = 



,2 yi2/2 



[r'+yinr^ + ymyl-yirdetgab. 



(2.34) 



Writing the Klein-Gordon field as 



itl'o'i/o „i4>l'^l „iJ/'2*i 



the Klein-Gordon equation {\/~gg^'^ dy^) — can be simplified to 



(2.35) 



1 



1 

^ ,,2\/,,2 



1 



I 1 /ai 



dr (r^-^ArdrRr 



ir^ + yim-yt) Ryi \yi 



Z5-6 



D-6 



(l-5'yf)(af-yf)(a^-2/f)9,,i?,, 



+ 



1 



+ 



(r2+y2)(y2_^2) j 

ir' + yl)ir' + y^) 

Ar 

{r^+yl){yl-yl) 



1 / \ 



ai 

^ ) (1 - g'yi){a{ - y22)(a^ - yi)dy,Ry, 

0-2 



D-6 



'(1) 

?3 



"(2) 
?3 



(r' + ylM-yl) 



"(1) 



j(j + i^-7) =0 



(2.36) 



where + D — 7) is the eigenvalue of the Laplacian on S ^. By putting in the values of i?^^,^ and by using the 
identities 



1 



1 



1 



1 



r2y2y2 (^2 ^y2)(-^2 ^y2)^2 (^2 _^ ^2) (^2 _ ^2)^2 (r^ + y|) (y2 _ y2 )y2 
^2 



2/i 



(r2+2;2)(r2+y2) ' (^2+y2)(y2_y2) ' (^2 + ^2) (^2 _ ^^2) 



1 



1 



1 



{r^ + yj){r^ + y^) {r^ + yf){yi-yf) (r^ + y^){yf - y^) 

we have the following separated equations 
1 



= 0, 



dr (r^-^ArdrRr 



Ryi \yi 



Rv2 \y2 



D-6 



■'VI 



D-6 



■'V2 



A 



1 / 4t 2t tn2 ajalj{j + D -7) 



{l-gVim-ytm-yt)dy,Ry, 



1 ^ 4t 2t t \2 ajalJij + D - 7) 
^(-2/f*o+2/?*i-«'2) -^2./u 



(l-g^y2')(a!-y2)(a2-y2)5,.i?,. 



-^2*0+ 2/2*1 -^'2) 



2 a 



?aij(j+I)-7) 



bir^ + 62 



biyj - 62 



(2.37) 
(2.38) 
(2.39) 

(2.40) 



(2.41) 



where bi and 62 are constants. 

In these equations the constants *i can be obtained by considering e*'^n*oe*'^i*ie*'''2*2 = e-^"*e™i'^ie*™2"^^ Using 
the relationship between t, 0i, 02 and t/iq: "01: ^2 in Eqs. (2.28) to (2.30), we have 



$2 



-u} + g'^{miai + 771202) , 

-uj{al + a\) + miai(l + .g^fla) -t- 771202(1 + .g^fli) 



-u)a\a\ + miaia\ + 77120^02 



(2.42) 
(2.43) 
(2.44) 
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FIG. 1: Due to the symmetries of the master equations we are able to patch together the QNM solutions obtained in the 
ai > a2 > octant to find the solution in the whole (ai, 02 )-parameter space. We are also able to find the angular eigenvalues 
61, 62 in this way. 



More explicitly, the radial equation is 



1 



pD-6 



dRr 

dr 



miOi 



m2a2 



- a? 



as 



and the angular equations are, for i — 1,2, 



777-202 



On 



alahU + D - 7) , 2 , 
^ ^ ^ , - bir^ - 62 



Rr = 



(2.45) 



D-6 



_d_ 



Z5-6 



{i-9M){<-yi){'4~yi) 



dRy 
dyi 



1 - (f-yl 

m\a\[\ - g^yf){al~ yf) _ ?T7^a^(l 

a{ - yf 



+ 2w[?T7iai(a2 - yf) + 771202(01 - yi)] - 2miOi7T7202(l - g'^y'f) 
'yf)iai-yf) alaUj + D-7) 



Vt 



yt 



biyf + 62 ^ = 



These master equations possess the following symmetries: 

(7711,01) O (7712,02) , and {m^,a.i) O (-771^,-0^) 



(2.46) 



(2.47) 



As we shall see, because of these symmetries, we only need to calculate the QNMs and eigenvalues in the Oi > 
02 > octant. The values in the other octants can be deduced from those with their quantum numbers transformed 
appropriately under the above symmetries, see figure 1. 

Assuming that the rotation parameters oi and 02 are small, the angular eigenvalues bi and 62 in the coupled angular 
equations above can be found perturbatively (i.e., as a power series in e = 02/01 and Oi), see Appendix A. The results 
are: 



bi 
b2 



Bi - 2uj{miai + 771202) + 2.g^ 771 1027^202 , 

2w(777iOi02 + Oi77l202) + 2r71iOi77l202 , 



alB2 



where 



Si = Bioo + Bw2 + O{e°ale^) , 
B2 = B200 + [S220 + B222 + Oie^t)] + 0{e^) 
The relevant terms can be obtained from equations (A15, A19, A27, A31, A32). 



(2.48) 
(2.49) 



(2.50) 
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III. NUMERICAL METHOD 



As shown in Appendix A a perturbative method can be used to determine the low order eigenvalues analytically; 
however, we were unable to do so in general for higher order terms, as exemplified in equation (A36). This is 
unfortunate because if approximate analytic expressions in terms of w existed for bi and &2 then we could have simply 
substituted them into the radial equation and performed the QNM analysis without any further reference to the yi 
equations (at least for small values of oi, 02). The fact that the perturbative expansions of bi and 62 are not algebraic 
expressions makes them unusable in the computation of the QNMs, because the analysis in general requires solving 
for the zeros of some polynomial equation in w. 

In addition to this, we can only expect the perturbative values to hold in the small rotation regime a2 < ai ^ 1. 
Therefore, it is desirable to have an alternative method of calculating these values. In this section we will describe a 
method that can be used to calculate both the eigenvalues bi, 62 and the QNM, w, numerically. This will serve as a 
consistency check of the results obtained in Appendix A and will also allow us to go to larger values in the rotation 
parameters. 

To achieve this we will use the improved Asymptotic Iteration Method (AIM) described in [31]. In the current 
problem this method has some advantages over that of the more commonly used Continued Fraction Method (CFM) 
[32]. In particular, the CFM requires lengthy calculations to prepare the recurrence relation coefficients that are 
subsequently used in the algorithm. Such manipulations are prone to error, and given the complexity of the current 
equations it is advantageous to use a method which bypasses this step. Furthermore, due to the existence of four 
Regular Singular Points (RSP) the CFM requires a further Gaussian elimination step in order to reduce the recurrence 
relation down to a three term recurrence relation [33]. As we shall see the AIM works for an Ordinary Differential 
Equation (ODE) with four RSPs in the same way as it would for a three RSP ODE and is therefore easier to 
implement.'^ 

The three equations (2.45) and (2.46) can be made to look more symmetrical by completing the square in terms of 
uj and defining: 



uj-{l + g^r^) 

^ - (1 - g'vl) 



7-2 + a\ 
miOi 



TO2a2 



TO2a2 



? 2 ' 2 2 / ' 



then, 



1 



j.D-6 



,D~6 , 



dr 



b,r' 



Rr 



(3.1) 
(3.2) 

(3.3) 



_d_ 
dyi 



D-6 



A, 



dRe. 



ial-yfnal-yf)\. , alaMj+D-7) 



Vt 



+ &i2/f-&2 Ue,(3.4) 



where i = 1,2. 

In the next section we shall solve the angular equations (3.4) showing how the AIM can be used to numerically find 
the bi and 62 angular eigenvalues and then we shall study the radial equation (3.3) and use the AIM to calculate the 
QNM, UJ. Before moving on, we shall briefly discuss the relation of ujr with super-radiance and the horizon structure. 



Super-radiance and the WKB form of the potential 



We can understand the form of cDj. when writing the radial equation in the WKB form by transforming as: 
where we defined the tortoise coordinate by 

dr^ _ (r^ + al){r'^ + al) 



dr 



(3.6) 



^ Nevertheless we found that the AIM required more iterations as the rotation parameter was increased which made it prohibitive to go 
beyond about ai ~ 1.5. Presumably the CFM would work more efficiently in this regime, however, in the current work only the small 
rotation QNMs were considered. 
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1 2 3 4 5 

FIG. 2: (Color Online) Plot of the D = 6 parameter space. The solid curve corresponds to solutions (in units of rn ~ 1) 
where the two horizons overlap (degenerate). Below the curve the outer horizon is fixed to unity. As the angular momenta 
are increased in the direction of the curve the two horizons cross (i.e., on the degenerate curve) and then pass into the blue 
shaded region corresponding to solutions with the inside horizon fixed to = I. Thus every solution below the curve (with 
normalisation such that the outer horizon is fixed at unity) has a corresponding equivalent solution above the curve with a 
different normalisation of the outside horizon. Due to symmetries we only need to study the 02 < ai region shown by the 
dashed line. Furthermore, in this work, we will only be investigating the small rotation ai < I region shown shaded in gray. 



The WKB wave equation is: 



d^Pr 

drl 

where is given in equation (3.1) and 

'alalj{j + D-7) 



ujI — 



(r2 + a\Y{r'^ + 



U{r) 



Uir) = 



r V(r2 + aj)^ (r^ + aj) 



2\2 



(l+.9V)(r^+af)(r^+a^) 3-Z?/2- 



+ as 



D-7 



2.9^ 



1 + g^r^ 



+ af 



A. 



3-D/2 1 



(3.7) 



"2' 
(3.8) 



Given the standard solution of the WKB wave function at infinity and near the horizon (mapped to minus infinity in 
tortoise coordinates), this identifies Cor{rh) < as super-radiant [34] for certain values of 01,02 and mi,m2. In the 
super-radiant case the transmission probability (|Ap) becomes negative and for Kerr-AdS we would expect this to 
lead to super-radiant instabilities, see [22] for an example in the simply rotating case. 

It may be worth mentioning that the WKB form of the radial potential could be used to find the QNMs via the 
WKB method of Iyer and Will [35] , once the angular eigenvalues are known [36] . It could also be used to find Hawking 
emissions via the WKB method [37] . 



Horizon structure for g = 

It is worth noting that even in fiat space (g — 0) the horizon structure is slightly different for D — 6 than it is for 
the larger dimensions. The reason for this is related to the number of angular momentum parameters that have been 
set to zero. In D = 6 there are only two possible parameters, yet in 13 = 7,8 there are three and in 13 = 9, 10 there 
are four parameters. In fixing the number of rotations to only two, as we have done in the current work, Z3 = 6 is the 
only case in which the full set are present. Furthermore, only D ~ 6 and D = 7 will have naked singularity solutions, 
see [38] for details. 

When obtaining our numerical results it will be necessary to fix a mass scale. Instead of setting M = 1 it is 
conventional (and convenient) to set the horizon radius = 1, see for example [20]. However, setting — 1 
automatically imposes the condition that a horizon exists. In this case, the solution will either have two horizons or 
a degenerate horizon. Using the results in [38], the degenerate horizon solution can be found in the even case from 
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the P2 polynomial. In units with = 1 the degenerate solution occurs when: 



A plot of this situation is shown in figure 2. As can be seen from this plot the degenerate horizon curve (solid line) 
divides the positive quadrant (oi, a2)-parameter space into two regions. One might wonder if the shaded blue region 
represents solutions with a naked singularity. However, this cannot be the case since the horizon was fixed to unity 
and therefore only those solutions with a horizon are being considered. The existence of the two regions is actually 
related to the fact that in general there are two horizons, an inner horizon and an outer horizon. It is the relative 
position of these horizons that is responsible for separating the parameter space into two. Even though we have set 

= 1 we have not specified which of the two horizons should take this value! Solutions under the degenerate curve in 
figure 2 correspond to those with the outer horizon fixed at unity while those above the curve correspond to solutions 
with the inner horizon fixed at unity. Since the two horizons are indistinguishable these two cases are really identical 
i.e., for every point above the line there is an equivalent solution below the line except that the normalisation of the 
mass is different in the two cases. Therefore, it is only necessary to study those solutions below the line in order to 
understand the entire stability problem. 

In a similar way to the above reasoning it is possible to show that there are no degenerate solutions (for only two 
non-zero spins) in higher than six even dimensions. Furthermore, using the Pi polynomial also defined in [38] one can 
again show that there are no degenerate solutions in odd dimensions*. Thus, Z? = 6 is in this sense a special case. 
Nevertheless, we will only be investigating the small rotation region (gray triangular region shown in figure 2), and 
in this regime the numerical method can be implemented identically in all dimensions. 



A. Higher dimensional spheroidal harmonics with two rotation parameters 

The two equations (3.4) are in fact the two-rotation generalisation of the higher dimensional spheroidal harmonics 
(HSHs) studied in [39]. In this case, the existence of two rotation parameters leads to a system of two coupled second 
order ODEs^. We note that, in general, one would expect that the generalisations of the HSHs to ['^^f-'^J rotation 
parameters would lead to even larger systems of equations. While these systems would also be useful generally in 
studies of MP black holes, here we will only focus on the two rotation case. 

It can be seen that these equations have regular singular points at yf = a\,a\,^ and 0. We assume that the 

cosmological constant is small and in particular that a\ <C \-^\- Recall that yi and 2/2 are defined on the domains 
shown in constraint (2.22). We would therefore expect the solutions to be well-behaved except possibly at the 
boundaries of these domains where singularities are present. In order to determine the regular solutions we need to 
define an appropriate norm on the space of solutions. First we change to the variable yf — The angular equations 
can then be written in the Sturm-Liouvillc form (assuming momentarily that uj and 62 are real): 

Au;(C,)^e.te) - - ^ (^^^^')^^^-(^»^) +9(e.)i?e,(6) (3.10) 
with the weight function wi(^i) = j^i^ ^ ^]-^g eigenvalue A = —bi, and 

Pte) = d'^^'^^'A?,, (3.11) 

where A^. and ui^^. are defined in the obvious way under the change of coordinates. Since w(^) > we can define the 



* This occurs even though in D = 7 there is a constraint on the angular momenta. In this case, the position of the degenerate horizon 
occurs at r/i = (this horizon has zero area and therefore should more properly be thought of as a naked singularity) and is therefore 
excluded by the assumption rj, = 1 since zero can not be scaled to 1 by an appropriate choice of units. 

^ For the moment we are considering u) to be an independent parameter. 
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two norm's: 

cx / 'er^'^^'l^eJ'dCi , (3.13) 



N2{Re,) (X / C2 " \Re.Vd^2 ■ (3.14) 







The rationale for this choice can be explained as follows. We can rewrite the Sturm-Liouville equation as an eigenvalue 
equation LR = XR where 



q{0 . (3.15) 



In analogy to the criterion discussed in [22] we note that for real lu and real 62, A (or bi) must be real. Therefore the 
inner product must be chosen so that L is self-adjoint when uj and 62 are real. From the Sturm-Liouville form it is 
easy to show^ that if the inner product is defined as: 

if, 9) - / riOgiOHOdi, (3.16) 

then L is self-adjoint i.e., {f,Lg) = {Lf,g). This inner product naturally induces the norms chosen above. 

However, one readily sees that the choice of norm is not unique. We could for example repeat the argument made 
above using A = 62 and in this case the weight function is found to be 1^2(^1) = The main point, however, 

is that even though the norms will give a different number (when acting on a given solution) , they will agree on which 
solutions are regular (finite norm) 

Under either choice of weight the regular solutions are found to be: 

i?i ~ (a-«^)^(a?-ei)^*i; ae(aia?), (3.17) 
^ e^/'(a2- 6)^*2; 6e(0,a2). (3_^8) 

Now for a given value of ui we can determine 61 and 62 simply by performing the improved AIM [31] on both of 
the angular equations separately. This will result in two equations in the two unknowns 61,62 which we can then 
solve using a numerical routine such as the built-in Mathematica functions NSolve or FindRoot. More specifically we 
rewrite equations (3.4) using (3.17) and (3.18) and transform them into the AIM form: 

d^^i , d^i ^ 

-TTT = Aoi— ^ + soi*i , 3.19 
d£,l d^i 

^ = Ao2^+S02*2. (3.20) 
"?2 "S2 

The AIM requires that a special point be taken about which the Agi and Sqi coefficients are expanded. As was shown 
in [40] different choices of this point can worsen or improve the speed of the convergence. In the absence of a clear 
selection criterion we simply choose this point conveniently in the middle of the domains: 

eor = ^, e„2 = #. (3.21) 



Eigenvalue results and comparison with double perturbation theory 



The numerical 61 and 62 eigenvalues for various parameters were computed and shown to be in good agreement 
with the perturbative values. Appendix A, for small e and lo. This serves as a consistency check between these two 



^ With appropriate boundary conditions. 

^ In Appendix A this ambiguity is somewhat more relevant. We find that in order to be able to simplify the expressions using the Jacobi 
orthonormality relations, one must choose the wi to normalize the Rg-^ solutions and W2 to normalize the Rg^ solutions respectively. 
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FIG. 3; (Color Online) D = 6, g — 0, [j,m\,m2,n\,n2) — (0,1,1,0,0). A plot of the eigenvalues for various choices of 
e = a2/a\. Note that the dependence on a\ has been scaled into the other quantities. 




FIG. 4: (Color Online) D = 6, e = a^jax = 1/2, (j, mi, m2,ni, 712) — (0, 1, 1, 0, 0). A plot of the eigenvalues for gai = 0.5i, 0, 0.5, 
corresponding to deSitter, flat, and anti-deSitter spacetimes respectively. Note that the dependence on ai has been scaled into 
the other quantities. 

methods. Some results are plotted in figures 3, 4. However, some issues arose and we found that both methods had 
their limitations, which we now briefly outline. 

We found that for e ^ 1 there was no appreciable difference between the numerical eigenvalues found after 16 or 
32 iterations. In other words the convergence was quite fast. As e — >■ 1 however we found that the convergence was 
much slower. For example, at ai = 3/2, 02 — 149/100 (i.e., e = 0.993) we needed about 80 iterations to get to the 
same level of accuracy that we required for smaller epsilon. See figure 5. In this case (for small values of oj) the 
perturbative method outperformed the AIM. 

However, we also found that the perturbative eigenvalues were very poor as lu became large. As an example, we 
choose the point ai = 1,02 = 1/2. Since e = 1/2 was relatively small we again found only 16 AIM iterations were 
required to get numerical convergence. However, with uj > 5 the perturbative values were clearly breaking down, see 
figure 6. 

The reason for this is that essentially the perturbation is expanded in the parameter cjoi where a; is assumed to be 
order unity. However, if oj is large then the error in this expansion becomes worse. 

The fact that the convergence of the numerical method is not sensitive to uj makes it more robust when calculating 
the QNMs, especially since we were only able to obtain the eigenvalues in the general case to 0{e'^,ai). Of course, 
what we have learned is that if we hold the number of AIM iterations at 16 then we would expect any values we 
calculate to have larger errors as e increases. In the next section we will calculate the QNMs completely numerically. 
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FIG. 5: (Color Online) D = 6, ai = 3/2, a2 = 149/100, (j, mi, m2, ni , 712) = (0, 1, 1, 0, 0) and g = 0. The dots are the numerical 
bi eigenvalues for an increasing number of AIM iterations, while the straight line is the value obtained from the perturbative 
method to 0{e^,al). 
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FIG. 6; (Color Online) _D = 6, ai = 1, 02 = 1/2, (j, mi, m2, ni, 712) = (0,1,1,0,0) and g = 0. The dotted blue line is the 
perturbation to order 0{e^,ai), while the dashed red line is the perturbation to order 0{e^,ai). The solid green line is the 
converged numerical bi eigenvalue. 

B. Radial quasi- normal modes 

Having successfully determined two independent methods for calculating the eigenvalues and having showed that 
they agree well with each other we can now proceed to calculate the QNMs with control over their range of validity. 
We start with the radial master equation (3.3). Thus far, when calculating the eigenvalues, we have been able to 
work in full generality, i.e., including the cosmological constant. However, due to the presence of new horizons and 
different boundary conditions, the flat, deSitter and anti-deSitter QNMs will need to be calculated separately. In this 
work we will focus only on the flat case, and from here on we set g — 0. 

Recall that QNMs are solutions to the radial master equation which satisfy the boundary condition that there are 
only waves ingoing at the black hole horizon and outgoing at asymptotic infinity. However, we found the AIM seems 
to work best on a compact domain. Therefore it is better to define the variable x — 1/r, so that infinity is mapped 
to zero and the outer horizon stays at = l/r^ — 1. The domain of x, therefore, will be [0,1]. Thus the QNM 
boundary condition is translated into the statement that the waves move leftward at x = and rightward at x = 1. 
We again choose the AIM point in the middle of the domain, i.e., at x = 1/2. 

In terms of x the radial equation (3.3) becomes: 

= -x--^l (--^-"^^f ) + ( ^"''^"^^If'-^"^^' ^^ - -l-ljij + D - 7)x^ -'^,-b.)R, (3.22) 

where Ax{x) = Ar(r = 1/x) and oJxix) = ujrii" = 1/a^)- 

After performing some asymptotic analysis, keeping in mind the definition (2.35), we find that for the solutions to 
satisfy the QNM boundary conditions we must have: 

i? - (1 - x)*"''"'-x(^-2)/2e*'^-/^y(x) , (3.23) 
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FIG. 7: (Color Online) D=6. Plots of the fundamental (j, mi, m2, ni, 712) = (0,0,0,0,0) QNM. On top are plots of the 
imaginary part and below are plots of the Real part. Left: a surface plot over the (ai, a2)-parameter space. Middle: a contour 
plot. Right: a plot of the QNM along the straight lines passing through the origin with gradient e shown in the graph. 



where 

= uj^ix^l) , (3.24) 
_ (l + a^)(l + al) 

= A^(x^l) • ^'-''^ 

We then substitute this ansatz into equation (3.22) and rewrite into the AIM form: 

y" - \oy' + Soy . (3.26) 

This final step was performed in Mathcmatica and then the resulting expressions for Aq and sq were fed into the AIM 
routine. 

The method we use to find the QNMs proceeds in a fashion similar to that used in [39, 41]. However, as already 
mentioned we use the AIM instead of the CFM. 

First we set the number of AIM iterations in both the eigenvalue and QNM calculations to sixteen*. We start 
with the Schwarzschild values (61,62,'^), i.e., at the point (01,02) ^ 0^ and then increment ai and 02 by some small 



Ideally this should be made as large as possible, however, we found that there was little difference in the computed QNMs when using 16 
or 32 iterations for ai < 1. This choice also gave good agreement in the small rotation regime with the perturbed eigenvalues calculated 
in appendix A. However, to go into the large rotation limit we found that much larger iterations were required to achieve convergence 
and this significantly slowed down the code. Thus the AIM method we have described here does not seem robust enough to explore the 
large rotation limits. 

Note we couldn't take the point (0,0) exactly as this would leave the yi and j/2 domains empty. To be precise we chose the point 
(ai,a2) = (1/50,1/100). 



15 




FIG. 8: (Color Online) D=6. Plots of the fundamental (j, mi, m2, ni, 712) = (0,1,0,0,0) QNM. On top are plots of the 
imaginary part and below are plots of the Real part. Left: a surface plot over the (ai, a2)-parameter space. Middle: a contour 
plot. Right: a plot of the QNM along the straight lines passing through the origin with gradient e shown in the graph. 



value"'^^. We take the initial eigenvalues (6i,&2)j insert them into the radial equation (3.23) then use the AIM to find 
the new QNM that is closest to to using the Mathematica routine FindRoot. 

We then take this new value of omega, w', insert it into the two angular equations (at the same value of ai and 
then solve using the AIM and searching closest to the previous 61 and 62 values. Thereby obtaining the new 
eigenvalues b'l, b^- We then repeat this process with the new (w', h'^^ 63) as the starting point until the results converge 
and we have achieved four decimal places of accuracy^^. When this occurs we increment ai and 02 again and repeat 
the process. In this way, we are able to find the QNMs and eigenvalues along lines passing approximately through 
the origin (i.e., starting from the near Scwharzschild values) in the (ai,a2) parameter space. We choose 6 straight 
lines with gradients of (e = 0, 0.2, 0.4, 0.6, 0.8, if^ and then use an interpolating function to interpolate the values in 
between these points. This then covers the ai > 02 > octant. 

Some preliminary results are shown for Z? = 6, in figures 7 and 8. We see that as we increase the gradient the 
imaginary part of the curves appear to be bounded between the 02 = to ai = 02 curves. Indeed, if this behaviour is 
a general phenomenon, then the most important regime for locating instabilities (i.e., when the imaginary part crosses 
the Im(uj) = axis) would appear to be the oi = 02 limit. Some general features of these plots are worth mentioning. 
For the case of vanishing angular modes mi_2 = 0, the solution is symmetric under horizontal and vertical reflections 



In our results we incremented ai by ^ at each step. The increment in 02 was then set by the gradient of the straight line taken in the 
{a\,a2)- parameter space 

We found that no more than 15 repetitions were required to achieve convergence. 

Note that if these lines went exactly through the origin, e = would mean a2 = and e = 1 would mean a2 = ai both of these situations 
would be pathological to the numerical method since the j/; domains would disappear. With the starting value at (1/50,1/100), we 
come very close to the single rotation and a\ = 02 cases for e = and e = 1 respectively while remaining in a valid domain of the 
numerical procedure. 



16 



Im [a>] 




Re [oj] 




1"^ -1 1 



FIG. 9: (Color Online) D=7. Plots of the fundamental (j, mi, m2,ni, 712) — (0, 1, 0, 0, 0) QNM. On top are plots of the imaginary 
part and below are plots of the Real part. Left: a surface plot over the (ai, a2)-parameter space. Middle: is a contour plot. 
Right: is a plot of the QNM along the straight lines passing through the origin with gradient e shown in the graph. 

in the ai and 02 axes in figure 7. Furthermore, with mi non-zero this reflection symmetry is broken in the real part 
of the QNM curves shown on the bottom right of figure 8, where they are skewed to the left. We have also confirmed 
this the case for m2 7^ 0. 

For higher dimensions such as Z? = 7, see figure 9, we observe similar behaviour to the D = 6 case where again 
we see the skewing of the real part for non-zero TOi , for example. The general dimensional dependence is shown in 
figure 10, where as is typical of singly rotating cases, larger dimensions lead to greater negative Im^ui) implying larger 
damping. These results also seem to indicate that larger values of D are more stable with increasing ai . 



IV. CONCLUDING REMARKS 



The intention of this work was to initiate the study of higher-dimensional Kerr-(A)dS black holes with with more 
than one rotation parameter for D > 5. In the present work we have considered those solutions with all rotation 
parameters set to zero except for two. As a first step in this direction, we have presented the general metric for such 
a spacetime. We have also separated the Klein-Gordon equation, writing out the corresponding radial and angular 
equations explicitly. In the general case with all rotations, e.g., see [8, 9, 42], the separation must be performed for 
each D separately. However, we found that in the doubly rotating case (ai 7^ 0, 02 7^ 0, 03 = 04 = • • • = 0) a general 
_D-dimensional expression could be obtained analogous to the one commonly used in the simply rotating case. 

It is worth stressing that in five dimensions there is only one one spheroidal equation, while in six and higher 
dimensions there are two angular equations, therefore in this work we only focussed on the D > 6 case (the five 
dimensional case will be considered elsewhere [30, 43]). We evaluated the QNM frequencies of the low- lying modes 
using a numerical AIM approach for both the angular and radial equations. In Appendix A, to get some quanti- 
tative understanding of the angular equations, we also developed perturbative expansions in powers of the rotation 
parameters e = a2/ai and ai for the angular eigenvalues. 
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FIG. 10: (Color Online) Dimensional dependence of fundamental (j, mi, m2, ni, 712) — (0,1,0,0,0) QNM for 02 = 0.4ai for 
< ai < 1. 



Our preliminary results for the QNMs suggest that slowly rotating black holes with two rotations are stable although 
our numerical code became slow for values of oi, 02 > 1, which unfortunately is also the region of most concern. More 
work in this direction, particularly for larger rotations and > (Kerr-AdS), would also be worthy of investigation. 

On the other hand, to discuss the stability of the ultra-spinning simply rotating black holes the angular eigenvalue 
in the large, imaginary, rotation limit [39], is typically relevant. For our case there are two rotation parameters. With 
one rotation parameter small and the other large, the situation will be very much like the simply rotating case [19, 20] 
and no instability is expected. Therefore it would more interesting to consider the other case with both rotation 
parameters large. This work will be pursued subsequently. 

In terms of other future work, with these separated equations we could also start to ask questions about the spectra 
of Hawking radiation (for five dimensions see [10]) and investigate the phenomenon of super-radiance in more detail. 
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Appendix A: Double Perturbation Theory 



In this appendix we develop expansions in powers of the parameters e = 02 /ai and ai for the angular separation 
constants bi and 62 defined in the equations (2.46). Since without loss of generality we have taken ai > 02, e < 1, 
furthermore, we assume Oi < 1. 

For brevity we shall outline the main steps and refer the reader to reference [30] for a more detailed account of the 
five-dimensional case^^. Since the latitude coordinates are restricted to a2 < yi < ai and < ^2 < 12, it is convenient 
to change variables to xi and X2 with 

yi = \ («i + ^1) -\ (a? - aj) xi ; 2/2 = \al (1 - X2) , (AI) 

where —1 < xi^X2 < 1. 

By making the convenient choice 

Bi = 61 + 2a;(miai -|- 771202) — 2(7^TOiaim2a2, (A-2) 

B2 = —2 [62 + 2w(r7iiaia2 + 0^771202) — 2miai?Ti2a2] (-^^3) 



For D = 5 there are two rotation parameters, but there is only one angular equation. This makes the pcrturbative analysis easier. 
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in Eq. (2.46) the perturbative expansion then only develops even powers of e and ai. In this case the operator 
equations are: 



and 



where 



[(1 - xi) + e2(l + X,)] {2 - g^al [(1 - x,) + 6^(1 + x,)] } 



\ {{D -5) + {D- l)xi) - \g'al{l ~ x,) {{D - 3) + p + l)xi) 



1 + a:i 



[{{D -5)-iD- l)xi) + {D + l)g^ajxiil - xi)] 

d 



2 VI -2^1 

eVa?(l + 2;i)\. . . 



alu^l + xi){l-e^) 



2\2 



4[2-g^al{{l-Xi) + e^l + x,))] 



~ xif + mle^jl + x^f] [2 - g^aj ((1 - xi) + 6^(1 + xi))] 
4(l-a;i)2(l + xi) 



(l-a:i)[(l-a;i)+e2(l + a;i)] 



O2 



4(1 -xi) 2(l-xi)j' 

i(l - xl) [2 - e\l - x,)] [2 - g^a\e\l - x,)] ^ 

^ [(i^ - 7) + (i^ - 3)X2] - ^(1 + - 2:2) [{D - 5) + [D - l)x2] 



rjg^alil - X2)^ [{D -3) + {D + l)x2] 

ilu;\^il + X2) [2-6^(1-3:2)] 
8[2-g^ale^l~X2)] 



dx2 



[2 - g^ajeHl - X2)] Tn\e\\ + X2f + (2 - £^(1 _ 3:2))' 



8(1 + 2:2) [2 -£2(1-^2) 
j(j + g-7) i?ie2(i-.x2) 
2(1-0:2) 8 



OiRi = -^Ri ■ 

We first expand the operators with respect to e. For each equation i = 1, 2 we have: 

(OiO + Ci2 + C'i4 + OiQ + • • • ) (i?io + Ri2 + ^'i4 + + ' ' ' ) 

= T (-BiO + Bi2 + -8^4 + Si6 + • • • ) [RiO + Ri2 + Ri4 + Ri6 + ' ' ' ) 



(A4) 



(A5) 



(A6) 



(AT) 



where the next subscript after the i refers to the power of e = 0-2/0-1 contained by those terms. Before going on we 
shall also need to look at the implications of the normalization condition to higher order terms. We first recall that 
the normalization conditions"'^* are given by: 



(l-xi)^i?? = 1, 



(1 -a;2)^i?: 



(A8) 
(A9) 



1'' See subsection (III A). 
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For convenience we define Ri = ^JWiRij2. Then schematically we have: 

y dx, [k^Q + Ra + ^,4 + + • • • ) (-Rio + + + + •■■)= 1 , (AlO) 
dxiii^Q = 1 ; / RioRi2 = ; /" dxi RioRt^ = — - / dxi i?^2 5 / ^2;^ RioRie = — [ dxi Ri2RiA ■ 



-1 



2 7-1 



1. Zeroth order in epsilon 0{e^) 
To the zeroth order in e, the eigenvalue equations are 



-^0 

OioRiQ = -^Rio, (All) 



where. 



^ ((i? - 5) + (i? - l)xi) - _ _ 3) + (i? + l)xi) 



A[2-g^al{l-x,)] 4(1 + 2(1 - Xi) J ' 

,2\ 



d 

dxi 
(A12) 



We first consider the ©20 equation. We note that the ©20 operator does not involve powers of oi at 0{e°) order, 
therefore we can write ©20 — O200 + C202 + • • • = C200 where the third subscript refers to the power of a, in those 
terms. To fix notation we also redefine the solution i?2o = R2Q0 and the eigenvalue B20 = i3200- 

One then observes that the corresponding eigenvalue equation 

^200 

0200^200 — ^ -R200 , (A14) 
is exactly solvable in terms of Jacobi polynomials. The solution is given by 

i?200(n2) = C2„,„,,,(l-X2r/2(l+X2)l"^l/2pif '^'''"^'^(X2), 

520o("2) = {2n2 + \m2\ + j){2n2 + \m2\ + j + D - 5) , (A15) 



where n2 = 0, 1, 2, ... . The normalization condition: 

(£>-7) ^ 

dx2(l-X2) 2 i?ioo(ni)i?ioo(ni) = (5„j„'^ , (AI6) 



1 



is satisfied if 
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_ / (2j + 2|m2| +4n2 + D- 5)r(n2 + l)r {j + 1^2! + + ^) | 

''"^"^^■"1 2i(2^+2|™.i+z.-7)r(|^,l + „, + i)r(, + n, + i^) / ■ ^""^'^ 

For the properties of Jacobi polynomials see [44]. 

Next we work on the Oiq equation. Following the same method we have used for e we expand the eigenvalue 
equation in powers of ai. The operator at zeroth order is 

/ o , d"^ 1 , _ , _ ^ yd ITT^'l , . N 



This can be found using the orthonormality of the Jacobi polynomials. 
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where, unlike O200, the Oioo operator is coupled to 0200 through -6200- 

The zeroth order eigenvalue equation C'iooi?ioo = | ^100^100 can again be solved and the solutions and eigenvalues 
are found to be: 



i?iooK) = Ci„,,„,(l-xi)^^(l + Xi)^Pi«^l"il)(a;i), 

Sioo("i) = (2{ni + n2) + \mi\ + \m2\ + j){2{ni + 712) + \mi\ + \m2\ + j + D ~ 3) , (A19) 
where ni = 0, 1, 2, ... , and for clarity we have defined 

= y S200 + {^^-J^) ' - 2n2 + |m2 \+j+l{D-5) . (A20) 

The orthonormality condition 

dxi{l-xi) 2 i?ioo(ni)i?ioo(ni) = (5„j„'^ , (A2f) 



4 

is satisfied with the normalisation 

_ / {2ni + K+\mi\ + l)r(ni + l)r {ni + k + |mi| + 1) ) 
'''"^"'^ 1 2«+l™i|-ir(ni + K+l)r(ni + |mi| + l) J ' ^ ' 

We now describe how given both the zeroth order eigenvalues and eigenfunctions we can go to higher order in the 
perturbative series. 

Second order in ai, 0{e , 

At next order in the perturbative expansion for Oi we find: 

C102P100 + Cioo^i02 = ^ (5io2-Rioo + -B100-R102) • (A23) 
Where the ©102 operator is found to be: 

O102 = lajg^l ~ - a;?)^ - Ja?5'(l - - 3 + p + + ^a?c^'(l + Xi) - ■ 

We also need to look at the effect of the normalization condition (A2f) to higher order terms. We again find a 
series of conditions analogous to (AfO), i.e., 

dxiRloo = 1 ; dxi i?iooi?io2 = • • • . (A24) 

-1 J-i 

Using these relations and the hermiticity of the operators we are able to find -B102 by contracting equation (A23) with 
^100: 

Bw2^ J dxi(l-xi)^i?ioo(ni)Oio2i?iooK) • (A25) 
Using the equation Oioo-Rioo = j^ioo^ioo to remove the second derivative we obtain: 



4 

r 1 . , 

C'102-RlOO = 



- ^a?.9'(l - 2;?)^ + ^alg^B2oo + + 2^1) - ^Biooa?<7'(l - a^i) i?ioo • (A26) 



After applying various identities for Jacobi polynomials [44] to remove the derivative and xi dependence in equation 
(Af 8) we obtain: 



j2 

S102 = -^{ij''+g'B,oo + 4.g'n,+g\2n + 5-D) + 2g'\mi\ 



(ml — k'^) 
(2ni + |mi| + K)(2ni + |mi| + K + 2) 



^(^2 ^ 2g^j^200 - g'B.oo + 9\2n + 5 - D) ~ 2g'\m,\) - 2alg'n,-^^-^^^^ . (A27) 
2 Zni + \mi \ + K 



21 



Note that in the single rotation limit (02 — >■ 0) with g = we find agreement with the result obtained in [39] using 
inverted continued fractions. The result does not, however, agree with that for g using inverted fractions in 
powers of Ci = aiu and ai = afg^ [40] because we used an expansion in different parameters: e = 02/01 and ai.^^ 

In theory wc could continue to go to higher order in a^, however, we stop at this order, to point out an issue that 
arises in the general case for the second order in e terms. 



2. Second order in epsilon ©(e^) 



We now consider the next order in e, firstly for O2 we have the operator equation, 

022^^20 + O20R22 ~ ^ (S20-R22 + B22R20) ■ 



(A28) 



After using the orthonormality of the doubly perturbed eigenfunctions up to second order in a\, the hermiticity of 
the operators and also the fact that many of the terms are simply zero in the O2 case, we find again that: 



S22i = / dX2{l-X2) 2 i?200C22i^200 , 



(A29) 



where i — 0,2. 

The operator takes the following form: 



O2 



-(1 + X2)(1-X2)2 {l + alg^) 



dx\ 4 



(1 - X2) (1 + alg^) {D-5 + iD- l)x2) 



dx2 



+ ^(1 + X2){aW + ttA) + ^^10(1 - X2) - (1 + alg^) 



2\ "^2(1 - ^2) 



4(1 + X2) 



Thus, working up to second order in a\, we have: 

' dx2 ' 8 



^220-^200 — 
0222-^200 = 



y (1 - a:^)^ + + X2) + ^Sioo(l - X2) + + ^^200(1 ~ ^2) 



R 



200 



- i(l - xDe^aig^^^ + ^-^uj\l + X2) + ^i?io2(l - ^2) + ^e'a?ff'j(j + ^ - 7) 



1 



e^a\g'^B2m{l - X2) 



i?2 



(A30) 



where in the above steps we used ©200^2001 cf. equation (A13), to remove second derivatives. Again using the 
functional properties of the Jacobi polynomials [44] we find: 



S220 = - (m^-Bioo + S2oo+4n2 + 2j + 2|m2|) 



(2^2 + |m2| + a)(2n2 + |m2| + a + 2) 



B222 — 



1 {ml - i?2oo + Bioo) + e^j{D + J - 7) - 2e^n2 , , ^"^f] 

2 2n2 + ITO2I + a 

i 

- (a?a;2 - B102 + alg^B2oo + 'ialg^n2 + 2alg'^j + 2alg^\m2\) 



+ 1^2!) , 



(A31) 



{ml - a^) 



(2n2 + ITO2I + a){2n2 + |to2| + a + 2) 
{a - [mal) 



^/'22,D 22d ^|222v7^|■ ■7\ o222 V" '"'2 ) ,222/- i |\ 

[a^uj + Bio2 - a^g B200) + e a^g j(D + j - 7) - 2e a^g n2- r- h e a^g {j - |m2|) , 



■2n,2 + |m2| + a 



(A32) 



where we defined a — j + {D — 7)/2. 



^® We have verified that perturbation theory (using orthogonal polynomials) for small ci, C2 and ai, 02 does indeed give the correct second 
order answer when C2,a2 — ^ for _D = 5 [30, 43] (cf. [40] using inverted continued fractions). 
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Second order in 0\ 

It is in the 0{e^) order of the Oi operator that our method runs into some difficulty in the general case. For the 
Oi operator at 0{e^,a\) order we find: 

Cl20^100 + C'l00-Rl20 ~ ^ (^100^120 + ^120-^100) i (A33) 

where after performing manipulations similar to those in the previous sections we find: 



/I 
(0-5) 
dxi{l-Xi) 2 i?iooOl20-RlOO 



(A34) 



As such, an expansion of ©12 in powers of a\ leads to 

n ^2 1 2.n ^ m 1^ .1^^ , e^mj (I + xi) 2 j(j+^-7) 



2(1 -xi) l-xi 4 

(ASS) 

Now we can use Oioo-Rioo, see equation (A18), to remove second derivative terms, and we obtain: 



(P-5) _ 



B120 = / rfxi(l-xi)^i?ioo -£^^-5)7^ 1^— +e 



(l + xi) d , 2(™2--B2oo)(l + a;i)-2j(j+i?-7) 



(1 — xi) dxi 2(1 — xi) 

B220 - e^ml 



2(1 



^100 ■ 

(A36) 



Unfortunately, the terms with factors of (1 — a;) in the denominator do not appear to allow any simplification via 
standard Jacobi identities [44]. Thus, in the general case it does not seem possible to find the solution to i?i2o in 
closed algebraic form. In principle one could still numerically integrate these expressions, however, as explained in 
section (III), when calculating the QNM's it is more advantageous to perform a full numerical calculation than to 
take this route. 

Nevertheless, we have found that the offending terms vanish for the special choice of parameters D — Q, j — {), 
mi = 7712 — 1 and 711=712= 0, and in this case, i?i2o ~ 0. 



3. Sixth order in the special case: D = &, j = Q, mi = = 1 and Tii = 712 = 

We also found that this special case allowed us to go to higher orders without encountering the above mentioned 
issue. To go to higher orders equation's like (A34) often can not be written simply in terms of the zeroth order 
eigenf unctions. In such cases we had to decompose the higher order functions in terms of linear superpositions of 
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zeroth order ones. We will give a more detailed explanation of this method in [30]. Here we just list the result: 



Bo 



20 

10 + ^ (2C.2 - 1752) al - — _ 53g2c^2 ^ 196^4)^4 



20 

H (2a; 

8444007^ 



- 645.9^0;^ + 28959/a;2 - 105644/)a? 

32 



^(2^2 - 17.g2)a? + ^(-' - 53gW + 196/)at - ^^^(-^^ - WT.g^c.^ + SlTSg^..^ - 184245^)a? 



2 + 6^ 



~8019^ 
20 

8444007 ^ 

2 + ^(40.2 _ 7^2)^2 - ^(0.4 - 53g2..2 _^ i9g^4)^4 



;(2a;'^ - 6455^^4 + 28959/0;^ - 105644/)a'J 



(A37) 



8444007 

(' 



- QAbg^uj"^ + 28959/^2 - 1056445^)a'j' 

~) - ^ ^ A \ A 32 /, 9 0\ / A 



8444007 



{2uj^ ~ 645g^uj^ + 28959/a;2 - 105644/)a? + • • • 



(A38) 



These results are compared to the AIM in figures 5 and 6. In figure 5 we see tl 
be as good as the perturbative value, where a similar plot was found for 62- In figure 6 we see 
eigenvalues break down for uj larger than ~ 5. 



and 6. In figure 5 we see that ~ 80 iterations are required to 

" — that the perturbative 



L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999), hep-ph/9905221. 
J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), hep-th/9711200. 
P. Kanti, Int. J. Mod. Phys. A19, 4899 (2004), hep-ph/0402168. 
R. C. Myers and M. J. Perry, Ann. Phys. 172, 304 (1986). 
R. Emparan and R. Myers, J. High Energy Phys. (2003). 
R. A. Konoplya and A. Zhidenko (2011), 1102.4014. 

S. W. Hawking, C. J. Hunter, and M. Taylor, Phys. Rev. D59, 064005 (1999), hep-th/9811056. 

G. W. Gibbons, H. Lu, D. N. Page, and C. N. Pope, Phys. Rev. Lett. 93, 171102 (2004), hep-th/0409155. 
W. Chen, H. Lu, and C. N. Pope, Class. Quant. Grav. 23, 5323 (2006), hep-th/0604125. 

H. Nomura, S. Yoshida, M. Tanabe, and K.-i. Maeda, Prog. Theor. Phys. 114, 707 (2005), hep-th/0502179. 
V. P. Frolov and D. Stojkovic, Phys. Rev. D67, 084004 (2003), gr-qc/0211055. 

M. Vasudevan and K. A. Stevens, Phys. Rev. D72, 124008 (2005), gr-qc/0507096. 

V. P. Frolov, P. Krtous, and D. Kubiznak, JHEP 02, 005 (2007), hep-th/0611245. 

P. Krtous, D. Kubiznak, D. N. Page, and V. P. Frolov, JHEP 02, 004 (2007), hep-th/0612029. 

V. P. Frolov and D. Kubiznak, Phys. Rev. Lett. 98, 011101 (2007), gr-qc/0605058. 

V. P. Frolov and D. Kubiznak, Class. Quant. Grav. 25, 154005 (2008), 0802.0322. 

A. Sergyeyev and P. Krtous, Phys. Rev. D77, 044033 (2008), 0711.4623. 

T. Oota and Y. Yasui, Phys. Lett. B659, 688 (2008), 0711.0078. 

Y. Morisawa and D. Ida, Phys. Rev. D71, 044022 (2005), gr-qc/0412070. 

V. Cardoso, G. Siopsis, and S. Yoshida, Phys. Rev. D71, 024019 (2005), hep-th/0412138. 

H. Kodama, R. A. Konoplya, and A. Zhidenko, Phys. Rev. D81, 044007 (2010), 0904.2154. 

H. Kodama, R. A. Konoplya, and A. Zhidenko, Phys. Rev. D79, 044003 (2009), 0812.0445. 

N. Uchikata, S. Yoshida, and T. Futamase, Phys. Rev. D80, 084020 (2009). 

V. Cardoso, O. J. C. Bias, and S. Yoshida, Phys. Rev. D74, 044008 (2006), hep-th/0607162. 

V. Cardoso, O. J. C. Bias, J. P. S. Lemos, and S. Yoshida, Phys. Rev. D70, 044039 (2004), hep-th/0404096. 

V. Cardoso and O. J. C. Bias, Phys. Rev. D70, 084011 (2004), hep-th/0405006. 

J. Boukas, H. T. Cho, A. S. CorneU, and W. Naylor, Phys. Rev. D80, 045021 (2009), 0906.1515. 

P. Kanti, H. Kodama, R. A. Konoplya, N. Pappas, and A. Zhidenko, Phys. Rev. D80, 084016 (2009), 0906.3845. 

P. Kanti, J. Phys. Conf. Ser. 189, 012020 (2009), 0903.2147. 

H. T. Cho, A. S. Cornell, J. Boukas, and W. Naylor, in preparation. 



24 



[31] H. T. Cho, A. S. Cornell, J. Doukas, and W. Naylor, Class. Quant. Grav. 27, 155004 (2010), 0912.2740. 

[32] E. W. Leaver, Phys. Rev. D41, 2986 (1990). 

[33] A. Zhidenko, Phys. Rev. D74, 064017 (2006), gr-qc/0607133. 

[34] B. S. DeWitt, Phys. Rept. 19, 295 (1975). 

[35] S. Iyer and C. M. Will, Phys. Rev. D35, 3621 (1987). 

[36] E. Seidel and S. Iyer, Phys. Rev. D41, 374 (1990). 

[37] A. S. Cornell, W. Naylor, and M. Sasaki, JHEP 02, 012 (2006), hep-th/0510009. 
[38] J. Doukas (2010), 1009.6118. 

[39] E. Berti, V. Cardoso, and M. Casals, Phys. Rev. D73, 024013 (2006), gr-qc/0511111. 

[40] H. T. Cho, A. S. Cornell, J. Doukas, and W. Naylor, Phys. Rev. D80, 064022 (2009), 0904.1867. 

[41] E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D 73, 064030 (2006). 

[42] G. W. Gibbons, H. Lu, D. N. Page, and C. N. Pope, J. Geom. Phys. 53, 49 (2005), hep-th/0404008. 

[43] A. N. Aliev and O. Dehce, Phys. Rev. D79, 024013 (2009), 0808.0280. 

[44] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic Press, UK, 2000). 



